Fluctuations and correlations in rotating Bose-Einstein condensates 
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We investigate the effects of correlations on the properties of the ground state of the rotating 
harmonically-trapped Bose gas by adding BogoHubov fluctuations to the mean-field ground state of 
an A^-particle single-vortex system. We demonstrate that the fluctuation-induced correlations lower 
the energy compared to that of the mean-field ground state, that the vortex core is pushed slightly 
away from the center of the trap, and that an unstable mode with negative energy (for rotations 
slower than a critical frequency) emerges in the energy spectrum, thus, pointing to a better state for 
slow rotation. We construct mean-field ground states of 0-, 1-, and 2-vortex states as a function of 
rotation rate and determine the critical frequencies for transitions between these states, as well as 
the critical frequency for appearance of a metastable state with an off-center vortex and its image 
vortex in the evanescent tail of the cloud. 



I. INTRODUCTION 

A rotating ultracold harmonically-trapped Bose gas is 
predicted to pass through many exotic phases with in- 
creasing rotation rate (for a recent review, see Ref. [1]). 
The mean-field description, which omits all correlations, 
predicts the zero-temperature ground state for a large 
number of particles to be a vortex lattice [21 13] , and is 
in good agreement with current experiments |3]. How- 
ever, exact diagonalization of the many-body Hamilto- 
nian [5] suggests the breakdown at very high rotation 
rates of the mean-field picture and a melting transition 
to strongly-correlated ground states, bosonic analogues 
of quantum Hall states [51 [7] . The onset of correlations 
and quantum fiuctuations can be expected to play a sig- 
nificant role in this transition to a strongly-correlated 
quantum liquid. However, a consistent theory of the zero- 
temperature melting of the vortex lattice does not exist 
so far [BHU] (for a theory of thermal melting of the lat- 
tice, see Ref. [I2])- A crucial first step in constructing 
such a theory is to understand better how correlations 
affect the system. 

With increasing rotation rate, the cloud expands in the 
transverse direction, and the particle density decreases. 
In each unit cell of the lattice, the vortex core occupies a 
larger fraction of the area of the cell [13] , and the aver- 
age displacement of the vortex from its equilibrium po- 
sition increases [^ due to the zero-point motion of the 
Tkachenko mode [HI [15] . Hence, the uncertainty in the 
position of vortices, which plays a leading role in the 
melting, increases at faster rotation rates [1]. This un- 
certainty is completely absent in the mean-field picture, 
in which the vortex positions are fixed and do not fluc- 
tuate. 

The nature of the correlations between particles 
changes as the rotation rate increases. For angular mo- 
mentum per particle less than or equal to unity (in units 
of h throughout), correlations in the exact ground state 
wave function jl6| are described by polynomials in the 
relative distance of the particles from the center-of-mass, 
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z ^ [x + iy) are the positions of the particles and Zc the 



center-of-mass in the complex plane. On the other hand, 
when the angular momentum per particle is of order the 
number of particles, correlations appear in the distances 
of particles from each other, as in the bosonic Laughlin 
wave function [T7] tp ^ Y\i<ij{zi — Zj)'^. 

The aim of the first part of this paper is to build rel- 
evant correlations on top of the mean-field many-body 
ground state and to investigate their effects on the ener- 
getics and physical properties of the ground state. Based 
on the inferred modifications of the ground state, the 
second part of this paper investigates, at the mean-field 
level, different ground states of the (0-, 1-, and 2-vortex) 
Bose gas and their respective transitions as the external 
rotation rate increases. 

Small-amplitude Bogoliubov fiuctuations around the 
mean-field ground state induce correlations by allowing 
small numbers of excitations to appear in nearby single- 
particle states. In a condensate with large number of 
vortices, the number of excited modes (single-particle 
harmonic oscillator eigenstates) involved is of order the 
number of vortices. Therefore, carrying out the general 
diagonalization is a mathematically challenging task for 
a many-vortex condensate. However, one can gain in- 
sight by working with a few-vortex system; for example, 
including the first three harmonic oscillator states is suf- 
ficient to describe systems with up to two vortices, as 
we show below. The simplest such system is a conden- 
sate with one singly-quantized vortex at the center of the 
trap, rotating at the critical frequency flc, at which the 
vortex becomes thermodynamically stable [HdS], and, 
hence, having unit angular momentum per particle. 

We find, indeed, that the correlations induced by Bo- 
goliubov fluctuations lead to a better ground state in the 
thermodynamic limit, lower in energy than the mean- 
field one. We also see that the fluctuations drive the 
vortex core away from the center of the trap by a fluc- 
tuating distance of 0{l/vN) (in units of the character- 
istic length of the trap). Moreover, for rotations slower 
than the critical frequency Qc, we find excitations with 
negative eigenenergy in the spectrum il8j which remove 
one unit of angular momentum from the gas, indicating 
an instability towards a lower-energy non-rotating state 
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FIG. 1: Schematic phase diagram of the condensate. 



for rotation rate n < flc and emphasizing the fact that 
the single-vortex mean-field ground state is not the best 
starting state. 

Based on these results, we construct a more energet- 
ically favorable mean-field condensate which, as a func- 
tion of ri, is either non-rotating, a single-vortex state, or 
a two-vortex state. The phase diagram in Fig. IT] summa- 
rizes our results. At a certain frequency flm (below flc), 
there exists a metastable state (a local minimum of the 
energy) with two off-center vortices in the cloud which 
are asymmetric about the origin; for a vortex close to 
the center of the trap, there exists an image vortex much 
further away where the particle density is negligible, in 
agreement with Ref. [T^. We calculate the critical fre- 
quencies, flc and fl2 respectively, at which the first and 
second vortices enter the cloud. The former agrees with 
the numerical result of Ref. [2] , while the latter is some- 
what larger owing to our incorporating only a restricted 
number of single-particle eigenstates in the mean-field 
ground state. 

In the next section, we outline the basic description 
of the rotating Bose gas in terms of Landau levels in 
the Coriolis force. In Sec. IIIIl we determine the small- 



amplitude fluctuations about the mean-field condensate 
and their effect on the properties of the ground state, 
and then in Sec. |IV[ we present the stable mean-field 
wave function that encapsulates the various phases of 
the rotating gas and their respective critical rotation fre- 
quencies. We set h = 1 throughout. 



II. BASIC MODEL 

We consider a cloud of TV bosons of mass m in a har- 
monic trap with frequencies w in the x-y plane and Uz 
in the z direction (with w^ ^ oj to tightly confine the 
gas in the axial direction), rotating around the z axis 
with angular velocity il. The characteristic oscillator 
lengths are d — l/y^mu in the transverse direction and 
dz = l/y^mujz in the axial direction. We assume weak 
two-body repulsive interactions of strength g — A-Ka/m 
where a is the s-wave scattering length. In the limit 
of fast rotation, fi < w, the gas becomes quasi-two- 
dimensional, and at zero temperature, it resides approx- 
imately in the ground state of the harmonic trap in the 
z direction. The many-body Hamiltonian in the rotat- 
ing frame is V.' = E^i(^° " ^^) + J2i<j 92D5{r, - r^), 
with the non-interacting single-particle Hamiltonian and 
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where r and p are the particle positions and momenta 
in the x~y plane, and 320 = g/v^dz is the effective 
interaction strength in two dimensions. 

The eigenstates of A" are the Landau levels \nm) - 
where the Landau level index n is the radial quantum 
number, and m > —n is the angular momentum along 
the rotation axis - with eigenvalues e^^ = (2n -|- |to| -|- 
1)0; — mVt. The characteristic energy scale of the two- 
body interaction is Vb = g2Til'2,TTd^ . The energy differ- 
ence between two successive higher Landau levels (n 7^ 0) 
is 0(2w) which, for fi — > w, is much larger than that 
0(0; — J7) between two states in a given Landau level. 
We assume the interactions to be sufficiently weak that 
Vb ^ 2w; thus, we ignore the small higher Landau level 
components in the wave function and safely assume that 
the system resides in the manifold of the n — Q lowest 
Landau level (LLL) states. With this assumption, the 
only relevant quantum number is the angular momen- 
tum index to; from now on, we drop the Landau level 
index n from the eigenfunctions, operators, and occupa- 
tion numbers for simplicity, unless otherwise noted. 

The LLL eigenfunctions and eigenenergies are 
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where z = {x + iy)/d is the position in the complex plane; 
we use r and z interchangeably. In terms of the corre- 
sponding creation and annihilation operators a]^ and am, 
the second-quantized Hamiltonian in the rotating frame 
is 

f" {rtii} 

where the interaction matrix element in the LLL is 
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III. BOGOLIUBOV HAMILTONIAN IN THE 
LOWEST LANDAU LEVEL 

We now turn to determining the effects of small- 
amplitude Bogoliubov fluctuations about the mean-field 
condensate on the properties of the system. We start 
with a condensate with a vortex at the center which ro- 
tates with angular frequency l^c, to be determined. We 
derive an effective LLL Hamiltonian along with its ex- 
citation spectrum (which includes an unstable normal 



mode) and show that its ground state has lower energy 
than the initial mean-field state in which all particles are 
condensed into the state |01). The initial condensate is 



4>{r) - VNiMr) 
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with iVi particles in |01), describing a vortex at the center 
with winding number 1. For A^i < N, we make the usual 
replacement of the operators a[ and Oi, corresponding 
to |01), by ^/]Vl in the hmit of large N. Although the 
total number of particles is fixed, interactions cause the 
number of particles in the condensate to fluctuate; the 
number of condensed particles can, thus, be written in 
terms of the total and the non-condensed particle num- 
bers as 



A^i =iV- 



.t 



(6) 



where the prime indicates that |01) is excluded from the 
sum. 

In the thermodynamic limit {N — > oo with iVVo con- 
stant), interaction terms that represent scattering of only 
one condensate particle or no condensate particles are re- 
spectively 0{1/\/N) and 0{1/N) smaller than those that 
involve two particles from the condensate and, thus, can 
be ignored. Following the standard procedure to write 
the Hamiltonian up to quadratic order in the excitation 
operators and using Eq. (|6| to conserve the total num- 
ber of particles, we derive the LLL Hamiltonian in the 
rotating frame, 

n' =[N{2uj -n) + iN^Vo] + [iNVo -{lo- n)]alao 
+ [iNVo + (a; - f7)]4a2 + -^NVo{alal + aoa2) . 

(7) 

The constant first term in square brackets is the energy of 
the mean-field state with all N particles condensed into 
|01) and no fluctuations present. The only two states 
in the LLL connected by the interactions in the pres- 
ence of a condensate in |01) are |00) and |02), i.e., for 
a LLL system, the maximum angular momentum trans- 
ferred in any scattering process is ±1, whereas allowing 
higher Landau levels brings in and connects |10) and |12), 
1 20) and |22), etc. Larger transfers of angular momen- 
tum take the system out of the LLL as well, e.g., |03) is 
connected to |1, —1) by a transfer of ±2 units, |2, —2) to 
1 04) by a transfer of ±3 units, etc. 

Conservation of angular momentum is reflected in the 
fact that any scattering process involves simultaneous 
transfers of +m and — m units of angular momentum 
(relative to the condensate) . The same method presented 
here was previously used by Linn and Fetter J18j to in- 
clude higher Landau levels perturbatively in the Bogoli- 
ubov excitation spectrum of this system. Also, Dodd et 
al. |20j have used a similar argument to describe angular 
momentum conservation in this system in the presence of 
an external perturbation. Later, Rokhsar reinterpreted 



their argument to show [3T] the existence of a negative 
energy excitation (the anomalous mode) with vortex core 
properties similar to those we find in Sec. |III A[ 

The canonical transformations to the bosonic quasi- 
particle operators 



a+i 



ua2 
uao 



V a 



(8) 



(with u and v real and positive) diagonalize the Hamil- 
tonian, provided that u^ = 2 and v'^ ~ 1. The new 
operators describe quasiparticles with ±1 units of angu- 
lar momentum relative to the condensate. Thus, the LLL 
Hamiltonian in the rotating frame becomes 

H' =[7V(2l^ - r!) + iTVVo - i^v^o] 

-|-(0-fic)alia-i + (w-0)a^ia+i, (9) 



where 



17, = w - ^ivyo- 



(10) 



This equation shows that the Bogoliubov ground state 
(with no excited quasiparticles) has lower energy com- 
pared to the mean-field one by — iVVb/4. Also, the nor- 
mal mode denoted by —1 has negative eigenenergy in the 
region f2 < J7c, indicating an instability in the system 
against being condensed into |01); this is the anomalous 
mode (see Ref . [22] and references therein) . Its existence 
shows that V'(r) is not the correct condensate for fi < fi, 
and, therefore, the Hamiltonian in Eq. ^k is not the cor- 
rect one for this regime. As the rotation rate increases 
beyond fie, further LLL states beyond {|00) , |01) , |02)} 
come into play in the ground state, especially once two 
or more vortices enter the cloud (see Sec. IV). Then, one 
must include Bogoliubov fluctuations around this new 
ground state in order to find the Bogoliubov Hamiltonian 
and its normal modes. For simplicity, we limit the dis- 
cussion here to = J7c which corresponds to our starting 
point, a system fully-condensed in |01). 

In the manifold of the first three lowest Landau levels, 
the field operator for removing a particle at position r is 



*(r) = (/)o(r) ao + 01 (r) a^ + 02(r) 02- 



Inverting Eq. ([8|) gives 



ao — ua-i — wck+i 
02 = wQf+i — V a_i 



(11) 



(12) 



and, thus, the expansion of ^(r) in terms of the quasi- 
particles is 



*(r) =7/;(r) + M</)2(r)a+i - W(?!)o(i")a+i 
-I- u0o(r)Q;_i - w02(r)a^i 



(13) 



This is the mode expansion »I'(r) = tp{'r)+^- [uj{r)aj — 
v*{r)aj] in terms of the quasiparticles (see, e.g.. 



Ref. [13]). Hence, the amplitudes for the —1 eigenmode 
are 



u_i(r) = u0o(r) 



(14) 



This mode, although having a negative eigenenergy for 
il < r^c, has a positive norm, since 



[\u-iirf - \v^i{r)\']d\ ^u'-v^ = l, (15) 



and is thus physical. Since (r |02)* — (r |2, —2), we arrive 
at the same amplitudes as derived up to zeroth order in 
Vq in Sec. Ill of Ref. ^8j. We note, however, that the 
only two states that are mixed, in fact, are |00) and |02) 
which are in the LLL, and not |2, —2) which is a higher 
Landau level; up to the level of the approximation used 
in this article, the fluctuations reside solely in the LLL. 



A. Properties of the Bogoliubov Ground State 

We now investigate the stable ground state of the 
Hamiltonian, Eq. ([9]), for fi = i7c and show that fluc- 
tuations drive the vortex away from the center of the 
trap and modify its velocity profile. The order param- 
eter is V'(r) = (^(r)) or, in terms of annihilation op- 
erators, (oi) — \fNi for the macroscopic condensate 
with A'^i particles in |01). The condensed state |iVi) is 
a coherent state that satisfies the eigenvalue equation 
fli \Ni) — y/Ni\Ni) with the normalized solution 



|iVi 



-JVi/2, 



/N^a] 



(16) 



where |vac) is the vacuum. This state docs not conserve 
particle number. 

The Bogoliubov ground state is determined by the con- 
dition that no quasiparticles be present, i.e., 

a±i |G) = 0. (17) 

Following the standard procedure (see, e.g., Ref. |24) ) . 
we find the normalized Bogoliubov ground state 

|G) = ^e-"^"«/^|A^i) (18) 

v2 

which has an expectation value of the angular momentum 
operator C — J2m''^^ln'^m given by (G| C |G) = A^. 

We now compare the lab- frame energy of |G) with that 
of the mean-field and exact ground states. The exact 
non-normalized many-body ground state for 2 < L < N 

is dg 

1p!;{zi ...Zn)= ^ {z^, - Zc){zi^ - Zc) • • • {Zi^ - Zc) 



il<i2<---<iL 



^N 



(19) 
where Zc = X^iLi ^i/^ is the center-of-mass coordinate; 
we suppress the factor exp[— ^^^j^ \zi~\ /2] common to 



all iV-particle LLL states from now on for brevity. This 
state has energy E^ = {N + L)oj + VoN{N - 1 - L/2)/2 
in the lab frame. The mean-field ground state for L = N 
with a vortex at origin. 



N 



Ipmfizi ...Zn)^ n0l(^j), 



(20) 



has energy E^^f = 2Nuj+V()N'^ /4 in the lab frame. There- 
fore, at L = TV (and, hence, at il — flc), the Bogoliubov 
ground state lies exactly half-way in energy between the 
mean-field ground state and the exact one. 

A diagnostic of the structure of the vortex is the cir- 
culation around a closed contour C encircling the center, 



r 



f v(r) 



dr, 



(21) 



quantized in units of 2Tih/m for a quantum vortex. 
The velocity is v(r) ~ (j(r)) / (/3(r)) where (j(r)) = 
(S/7TT.)Im<^^^^(r)V^(r)) is the expectation value of the 
current operator and (/o(r)) = ('l'^^(r)^(r)) is that of the 
density operator. In the mean-field state, Eq. (20), 



Vmt(r) 



h 9 

m r 



(22) 



which describes an irrotational superflow (except at the 
origin where the vortex is located) with circulation T^f — 
2Trh/m. However, the Bogoliubov ground state, (18 1, 
gives 



m 



N 1 

N 



(pW) = 



1 + iV 1- 



3 
2 3 



-Nl^ 



w' + ^ 



7rd2 

4-1 -\z\ 



e 



TT<P 



e, (23) 

(24) 



Thus, in the limit of large TV, the velocity field is 



VG(r) = — 



rO 



ni r^ + A(r) 



(25) 



where A(r) = (1 — \r'^)l(N + r^) is the correction due 
to quantum fluctuations. The circulation in |G) is then 



rG(r) = r,nf X y 



Nr^ 



r4 + Nr"^ + 1 



(26) 



For large but finite N , we find rc/Fnif ^ 1 (increasing 
from 1/2 to 4/3) for the large range \/^/N < r < \/]V 
while rG/r,„f — >■ as r — >■ and FG/Frnf — > 2 as r — > cx). 
The vortex (at the center of the trap in mean-field) is 
now pushed off-center by quantum fluctuations, hence 
the vanishing circulation as the contour shrinks towards 
the origin. The off-center vortex fluctuates very close 
about the origin as shown by the circulation approach- 
ing its mean-fleld value as the contour radius expands 
past 0{l/y/N). On the other hand, as the contour ex- 
pands even further towards infinity, the circulation grows 




|z| 



FIG. 2: Density (in units of l/ird^) of the Bogoliubov (solid 
line) and mean-field (dashed line) ground states for A'^ = 10, 
showing a vortex at the center. 



to twice its mean-field value, indicating the presence of 
an image vortex much further from the origin. These 
results agree with those of Sec. |IV[ fn the thermody- 
namic limit (TV — >■ cx)), however, Fq equals the quantum 
of circulation everywhere except at the origin (where it 
is zero) and at infinity (where it is twice the quantum 
of circulation); therefore, increasing number of particles 
suppresses quantum fluctuations of the vortex and leads 
to the vortex being driven less further from the center of 
the trap and the image vortex being driven more outward 
to infinity. 

Quantum fluctuations change the particle density of 
the ground state compared to its mean-field value, pmf — 
N \z\ e~l^l /ncP. The average density is now non-zero at 
the center of the trap, (p(0)) = l/nd'^ (see Fig. [2]). One 
can understand this finite density in terms of the vortex 
fluctuating about the origin, as discussed above. Snap- 
shots of the cloud in the laboratory would reveal a vor- 
tex at random positions (varying from shot to shot due 
to Bogoliubov fluctuations); averaging over these density 
snapshots would lead to Eq. (24 1 for the average parti- 
cle density of |G). One can also understand the finite 
density at the origin in terms of the single-particle quan- 
tum states, in particular, the non-zero occupation of |00) 
whose wave function does not vanish at the origin. 



Normalizing ip{z) leads to A^ = A^i 
and V are bounded by u^ + v"^ < N. 
rotating frame, then, becomes 



+ u'^ +v'^; hence, u 
The energy in the 



E' ^ 



N{2uj -n) + jN^Vq 



il — UJ -\ —— ]U 
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(28) 



The constant term is the energy of a condensate with N 
particles in |01). Note that E' is invariant under the si- 
multaneous transformation u — >• —u and v — >■ —v. We 
choose < V < u as this sector of the u-v plane is ener- 
getically favorable. 

We denote E' = E' - [N{2uj - Vt) + iV^Vb/4] as the 
energy contribution from the mixing of |00) and |02) with 
the condensate. Then, introducing the parametrization 
u — C,cosh{9/2) and v = Csinh(6'/2), we find 



J7 — uj_\_ + 



NV. 



NVa 



3 cosh 9 — v8 sinh ( 



E' ^e 

-f C^^ [-15 + 4cosh6l - 21 cosh(26') + 16\/2sinh(26i)] . 
128 

(29) 

Ignoring the quartic part for now, i.e., assuming u,w ^ 
v'^Vi or A^i < N ^ we minimize the quadratic part with 
respect to 9 and find tanh 9„i ~ a/8/3, which is depicted 
by the straight dashed line in the u-v plane in Fig. [3] 
With this value of 0, we have 



E' = {^~ r2c)C^ 



3T^ 
16 



C^ 



(30) 



The quadratic term shows that up to second order in 
the mixing due to interactions, the system is unstable for 
Q. <^c, as in the quantum treatment. 

For 51 > ilc, Eq. (30 1 is a monotonically increasing 
function of C with a minimum at C = or 



0, 



(31) 



IV. THE STABLE CONDENSATE 

The anomalous mode, denoted by —1 in Eq. ^, sug- 
gests that a condensate with a singly-quantized vortex at 
the center of the trap is not stable against fluctuations 
for fi < ric. However, Eqs. (13) and (14) indicate that a 



mean-field condensate wave function of the form 



V'(z) = y/N~i (j)i{z) + U (J)q{z) ~V(j)2{z). 



(27) 



can have lower energy and be stable, depending on the 
values of u and v (which can be taken to be real). We 
search for a better ground state by tuning the two extra 
degrees of freedom, u and w, as follows. 



describing a system fully condensed into |01) with one 
vortex at the center of the trap. (The superscripts "<" 
or ">" denote rotations slower or faster than U,^-) The 
energy in the rotating frame becomes 



E'> = (tvc^ + ^N^Vo) + N{n, - n) 



(32) 



where the first term is just the energy were all N particles 
condensed into |00). For 17 < Jlc, though, minimizing 
Eq. (pOl) gives 



'I6(17c-f^) 
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FIG. 3: Contour plot of E' for iV = 10^ iVK) = 0.1, and 
Q, = w — |A^Vb. Darker shades indicate lower energies. 
The straight dashed line represents the direction given by 
tanh^m = \/8/3 whereas the curved solid line is the solution 
of Eq. ([47|, discussed below. 



SO that u^ — ^/2v^, and the rotating-frame energy be- 
comes 



E„, — 



(Nu + -N^Vq) +N{nc -n)- 



4(f7c-rj)2 

5Vo 



(34) 



The boundedness of u and v implies < C™ cosh 9,n < 
A''. Hence, the region of validity of Eqs. (33 1 and (34 1 is 
flm < f2 < f2r where 



i Lr: 



UJ--NV0. 



(35) 



Then, for il < il^, the point (u<,u<) hes outside the 
circle defined by u^ + u^ = A^ and does not represent a 
physical solution. 

In order to check for the existence of lower-energy 
states on the edge of the circle, where u^ + v'^ = N, 
we have to compare Em with the corresponding energy 
Ee^ for points on the edge. Since iVi = on the edge, 
we find 



1 



3Ff 



0„.4 



E'^=[Noj+ ^N'Vaj + 2(r!, - ft^ + —v 



(36) 



For f2 < fie, this energy increases monotonically with v, 
with the minimum at 



u^ = VN, 



e = 0, 



(37) 



representing a system fully condensed into |00) (with no 
vortex) and the rotating-frame energy 



E'< :^Noj+^N^Vo. 
For fl > flc, however, the minimum is at 



< ^\ N- 



3% 



3Fn 



(38) 



(39) 



and the energy in the rotating frame is 



E'f = [Nuj + -N''Vq 



3Vn 



(40) 



Therefore, for $!„ < ^ < ^c, the difference between 
the two energies is 



EL 



E';. 



3Vr 



(nc-n){n-ni) 



(41) 



where Q,\ — uj — NV^. Since Q.\ < flm, we find that the 
point {u = \/N, V = 0) corresponds to the global mini- 
mum of the energy for the entire region < f2 < ^c, indi- 
cating that the system has fully condensed into |00) with 
no vortex, whereas a local minimum of the energy ap- 
pears at (u = u^, V = Vm) for flm < ^ < ^c- This latter 
point corresponds to a metastable state which describes 
two vortices (asymmetric with respect to the origin) in 
the condensate; hence, the metastability frequency O^ 
is the rotation frequency at which a metastable state ap- 
pears in the energy spectrum. On the other hand, for 
ri > fici the difference between the two energies is 



EL 



EL 



16 

31^ 



{n-nc){n-vQ 



(42) 



where Q.*2~uj~ j^NVq. Thus, iiVLc <^< fij, the point 
{u = Q,v — 0), i.e., the center of the circle, corresponds 
to the global minimum of the energy, indicating that the 
system has fully condensed into |01) with one vortex at 
the center. For fJ > 5I2, however, the global minimum 



of the energy is at {u 






on the edge of the 



circle, and the ground state is a coherent superposition 
of 1 00) and |02) with two vortices in a symmetric config- 
uration with respect to the origin. 

Hence, fie is the critical frequency for creating a cen- 
tered vortex, and U,2 is the critical frequency at which 
two vortices nucleate in the condensate. Note that fie 
agrees with the external rotation frequency derived in 
Refs. [2j[T8]; however, fij is bigger than the two- vortex 
nucleation frequency, w — 0.078A^Vb, calculated numeri- 
cally in Ref. [2] due to the very limited Hilbert space used 
here with only {|00) , |01) , |02)} as opposed to a rather 
large one used in Ref. [2]. 

The metastable state, for which u'^ — \/2'"^ 
vortices at the zeroes of ■(/'(•z), i.e., 



Vm, has two 



\m 



^N-i{v<Y±^N+{v<y 



(43) 



where $7 enters through v^ on the right side. Using 
Eq. ([33|, we find v< = \/Nji at Vl = f7„. Hence, 



\{Q,m) — ±v2. However, as 51 increases towards Jlc, 
we find for small v^ that z^ — > ±^/2(Vm/ ^/~Ny^^ or, in 
other words, z+ — )■ -|-(30 whereas z~ — >■ 0. As fi increases 
from J7m to J7c, the vortex at -|-v2 moves to infinity while 
the one at —y/2 reaches the center and becomes stable 
there. 



It is instructive to compare this result to that of 
Ref. [19^, where the authors assume a condensate with 
an off-center vortex at position b close to the center, i.e., 
x(^) ~ \N{z — b), and find perturbative corrections to 
the wave function using the Gross-Pitaevskii equation. 
The new non-normalized wave function in the LLL up to 
0(&2) is x{z) - VN{z -b){l + bz/2 + ■■■) (see Eq. (15) 
of Ref. HH]) with the rotation frequency, to lowest order 
in b, being fi ~ fie (see Eq. (16) of Ref. [H]). This wave 
function represents two vortices at b and —2/5. Normal- 
izing this wave function up to 0{b^) leads to 



X{z) 



-z 
2 



1 



\,^\= 



(44) 



We can repeat a similar procedure for the metastable 
state for which u.^j = \/2 v^ . Then, the condensate wave 



function, Eq. (27 1, becomes 



z + u,7, - 



< 



(45) 



Expanding for small u^ 
have 



and defining Ij — —u^J\^, we 



^{z) 



:z'+ l- 



(46) 



Thus, xiz) and ^^{z) have the same form with b -^ 6. 
Since 6 is small, ip{z) also represents two off-center vor- 
tices at S (close to the origin) and —2/S (much further 
away in the evanescent tail of the cloud). Using Eq. (33 1, 
we find the rotation rate of this two- vortex configuration 
in the lab frame to be fi = fl^ " (3/16) A^Vq \6\ which 
includes the next-order correction of 0(1^1 ) to the result 
of Ref. [TQ" . Note that while the calculations of Ref. [19] 
are limited to the fast rotating regime and are valid only 
in the vicinity of flc (due to their perturbative nature in 
the small parameter Q — uj), the method presented here 
covers the entire region < fl < uj and is only limited 
by the number of states included in the condensate wave 
function. 



which the change in the energy is extremum. We de- 
fine u = Rcosr] and v = Rsmrj and write 6E' = 
E'{u + Su,v + Sv) - E'{u,v) ~ {SuduE' + Svd^E'). 
We find the direction that extremizes the change in 
the energy by keeping R constant while varying 77; 
thus, Su ~ —Rsinrj 61] and 6v = Rcosr]5ri. Then, 
dE'/Srj = R{— sin rjduE' + cos rjdyE') = which gives 
dyE' /duE' = tan?7 = 6v/5u, where the last equality is 
just the slope of the tangent to the curve v{u). Therefore, 
the differential equation for the bottom of the valley is 



dv _ dyE' 
du duE' 



with 



v{u = 0) = 0. 



(47) 



Its solution ioT n — Lj ~ |iVVb is the solid line in Fig. k51 
We rewrite E' as C^^(^) + C^-B(^) where 



A{9) = ^^ [n + (3 cosh e-VS sinh 9)] 
8 

B(6i) = ^[-15-h4cosh0-21cosh(26')-|-16V2sinh(20)] 
128 



(48) 

26)] 
(49) 



with il = (O 



^NVo)/^NVo. At the critical points. 



d(^E' — deE' = 0. The trivial solution, C = 0, represents 
the center of the circle. Then, assuming C ^ 0, we find 



1 dA 
AlW 



1 dB 
2B~d0 



(50) 



which gives the critical points for all values of fi. 

For small J7, the only critical point in the valley is at 
the origin. However, Fig. |4] shows that as fl increases 
towards Qc, two other critical points (a saddle-point and 
a minimum of the energy) appear in the valley. It is also 
evident that the minimum moves towards the center as 
n increases. Thus, there should be a rotation frequency, 
namely the metastability frequency f2,„, at which the 
saddle-point and the minimum lie on top of each other 
and, hence, the second derivative of the energy vanishes. 
Using Eq. (50), we find that at ri„j. 



A. The Valley and the Metastable Point 

The derivation of the local minimum of the energy in 
the metastable regime, flm < ^ < ^c, has so far been 
restricted to small values of u and u, i.e., when J7 — > 
ric according to Eq. ( 33 ) . Ignoring the quartic term in 
Eq. (29) leads to a constant 0m, representing a straight 



line in the u-v plane which the local minimunr traverses 
as fl varies. Including the quartic terms causes the valley 
in the energy landscape to curve, as seen in Fig. [31 In this 
section, we rederive the metastable state and its onset 
frequency, Sim, for larger values of u and w, keeping the 
quartic terms in the energy. 

We first determine the equation governing the valley. 
The valley is a set of points, denoted here by v{u), at 



Id^A _ 1 d^B 
A1)(P ^ 2B1)(P 



1 dB 

2Bl)e 



(51) 



The single critical point (apart from the origin) at ^m 
satisfies both Eqs. ([50|) and ([5l[). Dividing Eq. ([5l]) by 



Eq. ( 50 ) leads to an equation for 6 independent of any 



other variable, namely 



d^A/dO^ d'^B/de'^ 



dA/de dB/de 



1 dB/de 

2 B ' 



(52) 



with the solution cosh6',„ ~ 2.0776 at Q,m or, using 
Eq. ( [50| ), ^m — —1.1824. Therefore, the frequency at 
which the first metastable state appears is, in fact, 



i in 



2-1824^^^^ 

: W NVq 



(53) 
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FIG. 4: Energy landscapes for fi = a; - (2.1824/8)iVVo, oj - (2.1094/8)iVK), and w - {l/4:)NVo from left to right, for iV = 10^ 
and NVo = 0.1. The dots represent the critical points with SP and MIN indicating the saddle-points and minima respectively. 




FIG. 5: The position of the two off-center vortices for N = 
10^, NVo = 0.1, and ilm < il < fie Note the different vertical 
scales for zi. and z~ . 



which is much closer to the critical frequency fie com- 
pared to the frequency given by Eq. ([35]). Contour plots 
of energy for rotation frequencies flm and flc can be seen 
in the left and right panels of Fig. |4J showing that the 
saddle-point at the origin (for fl < flc) turns into a min- 
imum for fl > ric- Since the approximations of the pre- 
vious section are valid for fi < fJc, the minimum ap- 



proaches the origin following the line v = u/v2, i.e., the 
slope of the valley at the origin is l/\/2 near fie- 

The metastable state located at (urmVm) represents 
two off-center vortices at the zeroes of the condensate 



wave function. 



where 



V2 



4i - V^ - ("m + %n) 4 



= 0. 



(54) 



The positions of the two roots of this equation as func- 
tions of n are plotted in Fig.[5]for the metastable regime, 
fljn < ^ < ^c- Just as before, one vortex approaches 
the center of the trap while the other moves to infin- 
ity as f2 increases. However, their initial positions are 
not symmetric with respect to the origin but are at 
z-{nm) = -0.5917 and 2+(f2„0 = -K4.2506 for the par- 
ticular values of N and Vq used in the figure. Since the 
vortices are stationary in the rotating frame, they pre- 
cess around the origin with frequency O in the lab frame. 
Therefore, as seen in the lab frame, z~ spirals in towards 
the center of the trap while 2+ spirals out to infinity as 
n increases. 
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